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ABSTRACT 

Using spatially-resolved spectra obtained with the Low Resolution Imaging Spectrometer at the 
Keck I telescope, we investigate the nature of ionizing sources and kinematic properties of emission- 
line gas in a LINER galaxy SDSS J091628.05-I-420818.7, which is a nearby {z = 0.0241) and bright 
{Mr = —20.2) early-type galaxy. After subtracting stellar absorption features using a combination of 
simple stellar population models, we measure the flux, line-of-sight velocity, and velocity dispersion of 
four emission lines, i.e.. Ha, H/3, [Oiii] A5007, and [Nil] A6584, to study radial change of emission-line 
fluxes and velocities. Compared to the point-spread-function of the observation, the emission-line 
region is slightly extended but comparable to the seeing size. The central concentration of emission- 
line gas suggests that ionization is triggered by a nuclear source, excluding old stellar population as 
ionizing sources. We find that emission-line gas is counter-rotating with respect to stellar component 
and that the [O ill] A5007 line is blueshifted compared to other emission lines, possibly due to an 
outflow. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: active — galaxies: kinematics and 
dynamics 



1. INTRODUCTION 

Active galactic nuclei (AGN) may play an impor- 
tant role in galaxy evolution by feeding energy to host 

filaxies, quenching star formation (e.g., iHopkins et al.l 
)05; Springcl ct al. 2005j) and leading to the black hole- 
galaxy scaJii]£_j;e]atiOT]sa^ 

verse (iFerra rese fc Me rrittI [2000t iGebhardt et al.l [2Q00l: 
iGiiltekin et al. 2009; Woo et al.l m\^ . By including 
simple AGN feedback prescriptions, recent semi-analytic 
models successfully reproduced several observed prop- 
erties includ ing luminos ity function, masses, and colors 
of galaxies (jCroton et a l. 2006; Bower et al. 2006), al- 
though it is still unclear how AGN feedback works in 
quenching star formation. Observational studies have 
found that the majority of early-type galaxies exhibit 
LINER (Low-Io nization Nuclear Emission-line Region, 
HeckmanI I1980D features at the cent ral region (e.g.. 
Ho et al.l [l997t lAnnibah etall l2010t ICapetti fc Baldil 



201H) . If LINERs are truly AGN, the high fraction of 



LINERs is consistent with the scenario that the AGN 
feedback works as a star-formation quenching mechanism 
in early- type galaxies. 

In spite of the large population of LINERs in the local 
universe, the nature of ionizing source of LINERs is still 
debated. On the one hand, LINERs have been gener- 
ally regarded a s a branch of low-luminos i ty AGNs (e.g. 
Ho et al.l Il997t iKauffmann et al.l 120031: iKewlev et atl 



20061) . This idea has been supported by several stud- 



ies on nearby LINER galaxies, which provided a num- 
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ber of smoking-gun evidences for AGN-origin of LINER 
features, e.g, the presence of broad emission-lines (e.g.. 
Hp ct al. 1997), compact emission-line region (e.g. iKeel 
1983: Pogge et al.. 20001 and u ltra-violet (UV) flux vari- 
abilitv fe.g.. lMaoz et al.ll2005D . 

On the other hand, narrow-line features of LINER (i.e.. 
Type 2 LINER) can be easily produced by other ion- 
izing mechanisms since the LINER features are char- 
acterized by strong low- ionization lines, such as [Oi], 
[Nil], and [Sil], relative to high- ionization and the 
recombination lines. Several studies suggested possi- 
ble mechanisms, e.g., evolved old stars, evolved young 
stars, and supcrnovac shock, for produci ng the LINER- 
like features (Bincttc ct al. 1994; Dopit a fc Sutherland! 
IT995I: iTa^uchi et al. 2000). By us ing ph otoionization 
models, for example, Binctte et al.l ()1994l) showed that 
evolved old stars, i.e., post-asymptotic giant branch 
(PAGB) stars, can produce Lyman continuum photons 
and exhibit LINER-like features, suggesting that a large 
fraction of LINERs among early-type galaxies can be 
origin ated by evolved stars (see also iStasinska et al.l 
I2OO8D . However, these models highly depend on the 
evolutionary tracks of PAGB stars, which are con- 
structed under a number of th eoretical assumptions (e.g., 
iStanghellini fc Renzinil [2000l ) . By comparing the mea- 
sured line flux ratios of a sample of nearby galaxies with 
hoto i onization mode l pred ictions from iBinette et al.l 
1994D . lAnnibah et al.l (|2010[ ) reported that 11 out of 49 
(~ 22%) LINER galaxies could be produced by PAGB 
stars, while the remaining ~ 78% galaxies required ad- 
ditional ionization m echanisms. S i nce th ey adopt the 
model prediction of IBinette et al.l (|1994l ). one should 
again k eep in mind the model uncertainties. In a simi- 
lar way, iTaniguchi et al.l ()2000D suggested young evolved 
stars, i.e., planetary nebulae, from recent starburst as 
a viable mechanism for LINER-like features, implying a 
close connection between post-starburst and LINER-like 
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phenomena. This scenario provides an ahernative ex- 
planation of LINER features particularly for late-type 
galaxies. In additoin to old and young stars, radia- 
tive shocks from oxygen-rich supernova can also produce 
LINER-like features without an AGN as suggested by 
iDopita k Sutherland (1.995). 

In contrast to the spatially concentrated LINER fea- 
tures originated from AGN ()Keell 119831 : iPogge et alj 
[2000t) . those LINER-like emissions due to non- 
AGN mechanisms are tend to be spatiall y ex- 
tended (fVeillcux ct al. 1995; Monrc al-Ibcro ct al. 120061 : 
Rich et airi2dl0l: ISaFzi et atl |2010|) . For example. 
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Pogge et al.( (|2000l ) investigated the size of narrow-line 



region (NLR) of 14 LINER galaxies based on the high 
spatial resolution narrow-band imaging data from the 
Hubble Space Telescope and showed that the most of the 
emission comes from central regions with the size of tens 
to hundreds of parsecs (pc) , indicating that central AGN 
activi ties are responsibl e for the LINER features. In con- 
trast, |Rich^F^ ()2010D reported extended LINER emis- 
sion in NGC 839, concluding that the LINER-like feature 
is related to other mechanisms than an AGN. 

It is also possible that LINER-like features observed 
with limited spatial resolutions may be a mixed bag o f 
various ionization mechanisms (e.g., iFilippenkol 119961 ). 
For nearby (z < 0.01) galaxies, spectroscopic observa- 
tions generally include nuclear emission from the cen- 
tral ^ 100 pc region, avoiding possible contamination 
from an extended emission-line region. In contrast, it is 
far more complicated to understand the nature of emis- 
sion lines in integrated spectra obtained with a large 
aperture; for exa mple, Sloan Digital Sky Survey (SDSS, 
I York et ani2000r P1 spectroscopic observations with a 3"- 
diameter fiber include emissions from much larger physi- 
cal scale (^ 5 kpc) in SDSS galaxies with a mean redshift 
z ~ 0.1. Hence, it is possible that even non-AGN galaxies 
can be interpreted as LINER galaxies due to contamina- 
tions from extended LINER-like emission region. While a 
number of studies classified SDSS LINERs as AGN (e.g., 
iKauffmann ct al. 2003; Kc wlev et al.ll2006[ ). other stud- 
ies argued that the S DSS LINERs may no t be related 
to AGN activity (e.g.. iStasinska et al.ll2008l: ISarzi et al.l 



.20111: lYan fc BlantonI l20Tli r 
2010r discussed that if SDSS 



[20T0I iCid Fernandes et al 
For example. ISarzi et al, 

LINER galaxies have modest equivalent width of [O ill] 
A5007 (< 2.4A), only a small fraction of these galaxies 
are truly power e d by AGN. A recent statistical study by 
lYan fc BlantonI (|2011[ ) further claimed that AGN is not a 
primary ionizing source in the SDSS LINER galaxies al- 
though they cannot rule out the presence of a true AGN 
hidden in the integrated spectra obtained with a large 
aperture in SDSS. 

Since the SDSS provides one of the largest galaxy sam- 
ple in the local universe, revealing the true nature of 
LINERs in the SDSS sample is a crucial step in under- 
standing AGN demography in the local universe, and 
the nature of the galaxy-AGN coevolution. Spatially- 
resolved spectroscopy can shed light on the nature of 
ionizing source by constraining whether the emission is 
from nuclear or extended region. The kinematic informa- 
tion of both stellar component and emission-line region 

^ |http: //www. sdss . org| 
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Fig. 1. — The AGN diagnostic diagram for the LINER early- 
type galaxies selected from the SDSS. Gray dots represent 1382 
nearby (0.02 < z < 0.05) and bright {Mr < -19.7) early-type 
galaxies with S/N > 3.0 for four emission lines, i.e., H/3, [Olll] 
A5007, [Ol] A6300, and Ho. A sample of 902 LINERs (black dots 
with error bars) are ind ctified by adopting an empirical demarca- 
tion line of lKewlev et al.i (2006.' ) between Seyferts and LINERs. We 
select a typical LINER early-type galaxies from the LINER sample 
(red point with error bar). The curved line indicates a theoreti- 
cal demarcation line between star-forming and AGN-host galaxies 
(Kewley et al. 20011). 

can also provide useful hints on the dynamical relation 
between emission-line gas and stellar component. 

To investigate the true nature of LINER features, 
we performed a pilot study using a carefully selected 
LINER early-type gal axy, SDSS J09162 8.05-I-420818.7, 
from SDSS DR7 (Aba zaiian et al.ll2009D . In this work, 
we present detailed analysis on the emission-line prop- 
erties based on the spatially-resolved spectroscopy. The 
paper is composed as follows. Section 2 describes sam- 
ple selection, observations and data reduction. Section 3 
presents analysis method on the extracted spectra. We 
present the main results on the spatial profile and kine- 
matics of emission lines in Section 4. Finally, we summa- 
rize and discuss our findings in Section 5. Throughout 
the paper, we use ACDM cosmological parameters, i.e.. 
Ho = 70 km s~i Mpc~\ r^A = 0.73, and = 0.27. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Sample Selection 

We selected a sample of low-redshift (0.02 < z < 0.05) 
and bright (r-band absolute magnitude Mr < —19.7) 
early-type galaxies with LINER features from SDSS 
DR7 to study the characteristics of LINER host galax- 
ies (Bae et al. in preparation). To represent the gen- 
eral population of LINERs which are mostly hosted by 
early-type galaxies and to avoid contribution from star- 
forming regions to the observed emission lines, we first 
selected early-type galaxies based on u — r color and 
Q — i color gradient using the criteria by IPark fc Choil 
(|2005D . Then, LINERs are identified following the cri- 
teria of lKewlev et al.l (I2006D on th e emission-line ratios 
(see Figure l: lBaldwiii etaLlflgsl . We used four emis- 
sion lines retrieved from the MPA-JHU catalog on the 
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Fig. 2. — SDSS optical (gri) composite color image of the selected 
LINER early-type galaxy, SDSS J091628.05-I-420818.7. The angle 
and the width of the black bar represent the position angle (-f 73°) 
and the width (1") of the slit, respectively. For comparison the 
3"-diameter SDSS fiber size is shown at the center (green circle). 

SDSS DR7 galaxiefl i.e., H/3, [Oiii] A5007, [01] A6300, 
and Ha, with signal-to-noise ratio (S/N) larger than 3.0. 
To examine whether recent star formation can account 
for LINER features, we further selected galaxies with the 
near-UV (NUV)— r color less than 5.5, which is an em- 
pirical demarcation line for the ea rly-type galaxi es with 
a recent (< 1 Gyr) star formation (jKavirai et al.| [2007). 

For a pilot study to constrain the nature of LINER 
emissions detected in SDSS galaxies, we chose from the 
sample one galaxy, SDSS J091628.05-I-420818.7, which 
represents general characteristics of the LINER early- 
type galaxies, with Mr = —20.2 and the effective radius 
of 2.7 kpc. At the redshift of 0.0241, the galaxy provides 
a reasonable spatial resolution, 0.51 kpc arcsec"^. Figure 
1 shows emission-line properties of the sample of SDSS 
LINER early-type galaxies with Mr < -19.7 and 0.02 < 
z < 0.05, as well as the selected galaxy. 

The galaxy's morphology in the SDSS optical im- 
ages indicates lack of recent and/or ongoing merger 
event (see Figure 2). Spatially extended UV-excess with 
NUV— r color was measured as 4.9 in the NUV images 
from GALEX (Galaxy Evolution Explorer. iMartin et al.l 
l2005rPI . suggesting that an existence of spatially extended 
young stellar population. No radio counterpart for the 
galaxy is det ected at 1.4 GHz through the VLA FIRST 
survey r,Becker et al.l [T995fl. 

2.2. Observations 

We performed a spatially-resolved spectroscopy on the 
selected target, SDSS J091628.05-I-420818.7, usi ng the 
Low Resolution Imaging S pectrometer (LRIS, Ok e et ahl 
119951 iRockosi et"aIll2010D at the Keck I telescope with 
1 hour (6 X 600 seconds) exposure on 31 March 2011. 
The D560 dichroic mirror was used to split the incoming 

^ |http : //w ww ■ mpa-garc hing ■ mpg . de/SDSS/DR7 1 
^ http: //galex. stsci ■ 6du| 
* http: //sundog. stsci . edu] 



light into the red and blue side CCDs, and the 600/4000 
grism and the 600/7500 grating were used for the blue 
and red sides, respectively. The instrumental setup was 
chosen to cover the spectral range from 3200A to 9000A, 
including forbidden hues [Oiii] A5007, [Nil] A6584, and 
Balmer lines, i.e.. Ha, H/3, which are essential ingredients 
to identify LINERs. We used a l"-wide slit with the 
position angle of -1-73° counter-clockwise direction from 
the North to the East, which was determined as a major 
axis based on /i?^F/ellipse fitting using the SDSS g- 
band image. To increase S/N, we used 2x2 binning in 
the CCD, resulting in a spatial scale 0."27 pixel"^. The 
dispersion scales for the blue and red sides are 1.23A 
pixel^^ and I.6OA pixel"^, respectively. 

To calibrate the flux and estimate the seeing size, 
we obser ved two sp ectrophotometric stars, Feige 34 and 
Feige 66 (lOkelll99CT , right before and after observing the 
target galaxy. The seeing size decreased during our ob- 
servations from ^ l."6 to ~ l."l and from ~ l."2 to 
^ 0."9 on the blue and red side, respectively. The effect 
of seeing and its variation is discussed in Section 3.3. 

2.3. Data reduction 

We reduced the spectroscopic data using a series of 
I RAF scripts inc luding ^wodspec p ackage as pre viously 
performed (e.g., IWoo et al.l [20061: IMcGill et all 12008.) . 
The reduction process includes standard procedures, i.e., 
bias subtraction, flat-fielding, wavelength and flux cali- 
brations. After bias subtraction and flat-flelding, we re- 
moved cosmic rays u sing the Laplacian cosmic ray iden- 
tiflcation algorithm (jvan Dokkum 200 L). Arc lines (Hg, 
Ne, Ar, Cd, and Zn) and sky emission lines were used to 
calibrate wavelength scales for blue and red sides, respec- 
tively. Atmospheric extinction was corrected for each ex- 
posure while the galactic extinction was corrected using 
deredden task. Flux-calibration was performed using two 
spectrophotometric standard stars Feige 34 and Feige 66, 
which were observed with the same instrumental setup 
during the observing run. After subtracting sky back- 
ground by adopting sky regions at > 20" from the galaxy 
center, we extracted nine spectra consecutively along the 
spatial direction using a 4-pixel aperture diameter, which 
is equivalent to l."08 and ~ 0.55 kpc at the distance of 
the target. However, we only used flve central apertures 
in our analysis, since we could not detect weak emission 
lines, especially H/3, from the spectra extracted from the 
outer region. Thus, we cover the spatial range from —1.1 
kpc (SW) to -t-1.1 kpc (NE) with respect to the center 
of the galaxy. 

Finally, we combined the extracted spectra from each 
exposure with a median filtering, after scaling fiuxes to 
a mean flux averaged over from 4950 A to 5 150 A and 
from 6400 A to 6800 A for blue and red sides, respectively. 
The selected ranges are most preferred since they include 
useful emission lines for AGN diagnostics, e.g., H/3 and 
[Oiii] A5007 on the blue side, [01] A6300, Ha, and [Nil] 
A6584 on the red side. Using arc and sky emission lines, 
we measured the instrumental resolutions as 6.1A and 
4. 2 A full width at half maximum (FWHM), respectively 
for blue and red sides, corresponding to Gaussian veloc- 
ities ~ 152 and ~ 80 km at the wavelength of the 
observed H/3 and Ha. The final reduced spectra along 
the spatial direction are presented in Figure 3. 
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Fig. 3. — The observed spectra at each apoture along the sUt (black line) are plotted over the best-fit stellar models (red line) at each 
panel. The corresponding distance from the center is denoted in the upper-left Conner. The spectral range of each emission line, i.e., H/3, 
[Olll] A5007, and Ha-(-[Nll] doublet, is makred with green shades from shorter to longer wavelength, respectively. All spectra are shifted 
to the rest frame using the velocity measured from the central aperture. 



3. ANALYSIS 

In this section, we describe the method for measuring 
the emission-hne fluxes. Since it is crucial to remove stel- 
lar absorption features to reliably measure the emission 
line fluxes, we first present the procedure for subtracting 
stellar absorption line (e.g.. Ho 2008). Then, we present 
the emission line flux measurements. Lastly, we describe 
how we construct the point-spread- function (PSF) . 

3.1. Stellar Continuum Subtraction 

In order to subtract the stellar absorption lines, we 
need to construct a best-fit stellar model with a combi- 
nation of stellar templates or population models. Here we 
exploit 47 MILES simple stellar pop ulation (SSP) tem- 
plates (jSanchez-Blazguez et al.lf2OO60 . which include the 
population age from 60 Myrs to 12.6 Gyrs, with twice- 
solar met allicity (Z = 0.03) and the Kroupa initial mass 
function (|Kroupall2001l) . For a consistency check we com- 
pare the results using simple stellar population models 
with solar metallicity {Z — 0.02), leading to little differ- 
ence with slightly higher value of the fitting. 



We construct a best-fit stellar model for galaxy 
spectra using the penalized pixel-fitting code (pPXF, 
ICappellari fc EmsellemI [2004r Fl. which is widely used in 
stellar continuum subtraction and gas-kinematics anal- 
ysis for individual galaxie s and large surveys (e.g.. 
iFalcon-Barroso et al . 2004; Cappellari et al.ll2011l). Th e 
pPXF implemented the MPFIT code (Markwardt 2009), 
which is ba sed on the Levenberg-Marquardt algorithm 
(jMorel[r978| ) for minimization. For the fitting pro- 
cess, we use the limited rest-frame wavelength range from 
3550A to 7380A, to match the available wavelength range 
of the MILES SSP templates 

We mask out the wavelength ranges of well-known neb- 
ular emission lines and sky absorption features during the 
fitting process. To account for imperfect flux calibration 
and increase the quality of the fit, we adopt a multiplica- 
tive polynomial function. Note that we choose a high 
order (100-200) polynomial depending on the shape of 
spectra, however the polynomial function does not sig- 

^ http:/ /www-astro. physics. ox. ac.uk/~mxc/idl 
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nificantly affect the results since the polynomial func- 
tion does not attempt to fit any single emission line. 
During this process, we simultaneously measure the first 
(line-of-sight velocity) and second moments (velocity dis- 
persion) of the stellar absorption lines using spectra ob- 
tained at each aperture, which are used to construct the 
radial profile of rotational velocity and velocity disper- 
sion of stellar component as described in Section 4.2. To 
measure stellar velocity dispersion, we use the '-^ 200A 
range of spectra including the Mg& triplet and Fe (5270, 
5300A) lines. We correct for the instrumental resolutions 
of the observed spectra and the MILES SSP templates 
(FWHM=2.5lA) by subtracting them in quadrature (see 
Eq. 1 in Woo et al. 2004). The typical measurement er- 
rors of stellar velocity dispersion determined from the fit- 
ting procedure are ~ 5%. Using the spectrum extracted 
from the centeral aperture, we measure stellar velosity 
dispersion as 178 ± 10 km s~^, which is consistent with 
182 ± 3.8 km sec~^ measured from the 3" diameter fiber 
spectrum of SDSS. In Figure 3, we present the best-fit 
stellar models overlaid with the observed spectra at each 
aperture. All spectra are shifted from the observed-frame 
to the rest-frame using the radial velocity measured from 
the central aperture. 

After the continuum fit, we subtract the stellar model 
from the observed spectra, revealing the stellar absorp- 
tion corrected nebular emission lines, which are used 
in emission- line measurement (see Section 3.2 and Fig- 
ure 4). The standard deviations of the residuals from 
stellar subtraction (excluding emission-line regions) are 
7.0 X 10~^^ erg cm"^ s~^ A"^ and 1.8 x 10"^* erg cm"^ 
s""'^ A~^ for apertures at the center and outer (-1-1.1 kpc), 
respectively. For comparison, peaks of flux density for 
each emission line are significantly larger than the resid- 
ual from stellar subtraction (> 3(t), confirming that the 
stellar fit is good enough to measure emission-line fluxes. 
In case of the 11/3 at -1-1.1 kpc aperture, which is the 
weakest emission-line we measured, the peak of fiux den- 
sity is comparable to the fitting residual, indicating that 
it is difficult to measure the weak H/3 emission-line fiux 
at -fl.l kpc aperture. 

3.2. Emission-line Measurement 

After the stellar subtraction, we measure four emission 
hues for each aperture; II/3, [Oiii] A5007, Ha, and [Nil] 
A6584. The wavelength range we used also contains the 
other emission lines, such as [S ii], [O i], that are typically 
used in AGN diagnostics. However, we could not use 
those lines since the [Sii] AA 6716, 6731 hues fall in the 
atmospheric absorption B band and the [O i] A6300 line 
is too weak to properly measure line flux. In Figure 
4, we present the emission lines after subtracting stellar 
component and the best-fit single Gaussian models. 

We use the Monte-Carlo technique to measure the 
mean strength, wavelength, dispersion and estimate their 
uncertainties for each emission line simultaneously. First 
of all, we take a standard deviation of the continuum fit- 
ting residuals within ^ ±100A range from the center of 
each emission line as a proxy for systematic uncertain- 
ties, which is larger than the measurement uncertainties 
in the emission-line fitting process. Thus the uncertain- 
ties adopted in the following simulations should be con- 
sidered as upper limits. Then, we generate a thousand 
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Fig. 4. — Residual emission lines (gray line) after subtracting 
stellar continuum and the best-fit emission line model (red line) for 
the H/3-f [Oiii] (left panel) and the Hq:+[Nii] regions (right panel), 
repectively. Vertical lines represent the reference wavelength of 
each emission line. The distance from the center is denoted in left 
panels nd side. Orange lines on the right-hand side panels exhibit 
the sum of three Gaussian functions for Ha and [N ll] doublet. 

simulated spectra using the random values at each pixel, 
which follow a normal distribution with the systematic 
uncertainty as a standard deviation, and add them to the 
emission lines. For each simulated spectrum, we measure 
line fluxes using a Gaussian model fitting method as fol- 
lows. For the Ha-|-[Nii] region, we use triple Gaussian 
models by constraining the [Nil] A6584-to-[Nii] A6548 ra- 
tio as 3.0, while we use a single Gaussian model for other 
emission lines. In this procedure, we measure the central 
wavelength and line dispersion and their uncertainties 
of each emission line simultaneously, in order to investi- 
gate the kinematics of nebular emission. We correct for 
the instrumental broadening by subtracting the instru- 
mental resolution from the measured line dispersion in 
quadrature as ge nerally practiced (e.g.. lBarth et al.|[2002l : 
iWoo et al.|[200l . In case of the [Oiii] A5007 at 1.1 kpc 
aperture, the measured line dispersion (~ 2.4A) is com- 
parable to the instrumental resolution 2.5A for the 
blue side), hence we did not correct for the instrumen- 
tal resolution and the measured line dispersion of [O ill] 
should be taken as an upper limit. 

Systematic uncertainty of the radial velocity comes 
from the uncertainty of wavelength calibration, which 
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Fig. 5. — PSF models based on the observation of Feige 66 for 
H/3 (blue side; left panel), and Ha (red side; right panel). The 
photon counts measured at each pixel along the slit are normalized 
by the peak value (blue dots). Two Gaussin functions (red lines) 
are used for the best-fit PSF model (black line). The residuals are 
shown at the bottom of each panel (dashed line). 

was estimated as~ 0.2 A by adopting the standard de- 
viation of arc or sky line wavelengths with respect to 
the wavelength solution. The measurement errors of the 
line velocity derived from the Monte-Carlo realization are 
less than or comparable to the systematic error in most 
cases. Hence, we added both errors in quadrature to ob- 
tain the final uncertainty of each emission-line velocity. 
For the line dispersion measurements, on the other hand, 
we only consider the measurement error derived from the 
Monte-Carlo realization. 

3.3. Point Spread Function 

A precise measurement of a PSF is an important task 
to examine whether the LINER-like features are origi- 
nated from an extended or nuclear region. To measure 
the PSF we utilize the spectrophotometric stars observed 
right before and after observing the target galaxy. Using 
photon counts of each pixel from standard star observa- 
tions, we obtained the spatial profile of PSF at each wave- 
length range which is corresponding to the spectral range 
of each emission lines. For each emission-line region, we 
performed a median combine by using 10 pixels along the 
dispersion direction, then we constructed a model PSF 
by fitting the observed PSF with multi-Gaussian func- 
tions, resulting in one PSF for each emission line. 

Figure 5 illustrates the PSF model based on the spa- 
tial distribution of a standard star Feige 66. Note that 
the spatial distribution is slightly extended than a sin- 
gle Gaussian. Thus, we used two Gaussian components 
for the best PSF model. In the H/3 region (Fig. 5 left), 
we used FWHM 4. "91 and l."08, respectively for weak 
and strong Gaussian compoents with the peak ratio of 
weak to strong components as 0.03 while in the Ha region 
(Fig. 5 right) FWHM was 4. "91 and l."08, respectivley 
for weak and strong Gaussian components with the peak 
ratio of 0.06. 

4. RESULTS 

In this section, we present the results of spatial profiles 
of emission-line flux and kinematic analysis for emission- 
line gas and stellar component, and suggest the most 
probable ionizing sources for the observed LINER fea- 
tures, and possible explanations for the observed kine- 
matics. 
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Fig. 6. — Spatial profiles of emission-line flux: [Olll] A5007 (top 
left), H/3 (bottom left), [Nil] A6584 (top right), and H« (bot- 
tom right). Thick solid lines indicate the profile of each emission 
line, and vertical error bars indicate a standard deviation of the 
measured fluxes from the Monte-Carlo realization, while horizon- 
tal error bars represent the size of aperture for each emission-line 
flux. Dotted lines are normalized stellar continuum fiux at 5100A 
(F5100). Gray shade for each panel indicates that the minimum- to- 
maximum PSF range of the observation. Wc also plot the mean of 
both PSFs along the slit with a thin solid line over the gray shade. 

4.1. Spatial Profiles of Emission-line Flux 

To examine whether the emission-line region is spa- 
tially extended, we investigate the radial distribution of 
each emission line by comparing the spatial profiles of 
emission lines with that of stellar component and the 
PSF determined in Section 3.3. Figure 6 presents the 
spatial profiles of the emission-line flux for [O ill] A5007, 
H/3, [Nil] A6584, and Ha along the slit. We overlay the 
PSF determined in Section 3.3, after normalizing it to 
the emission-line flux at the center. Since the exact PSF 
at the time of exposure on the galaxy is not known, we 
assume the true PSF size is between the maximum and 
minimum PSF sizes of spectrophotometric stars as de- 
termined in Section 3.3, as plotted with black solid lines 
and gray shades. 

We find that the spatial profiles of emission-line flux 
are slightly more extended than, but comparable to the 
PSFs, indicating that the emission-line gas is centrally 
concentrated. At —1.1 kpc, there is a noticeable excess 
of emission-line fluxes compared to the PSF, suggesting 
that the emission-line region is slightly more extended 
than a point source. We compare the spatial profiles of 
emission- line region and stellar component, to investi- 
gate a physical link between stars and ionizing sources of 
the emission lines. In Figure 6, the spatial profile of stel- 
lar continuum flux at 5100A (i^sioo) is presented after 
normalizing it to each emission-line flux at the center. 
Compared to the spatial distribution of stellar compo- 
nent, the emission-line region is centrally concentrated, 
implying that old stellar population may not be a pri- 
mary ionizing source of the LINER feature. Unless old 
stellar population is concentrated at the galaxy center. 
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Fig. 7. — The AGN diagnostics based on the emission-line ra- 
tios for each aperture. Compared to the emission-line ratio ob- 
tained from the SDSS integrated spectrum (red dot), the emission- 
line ratios measured from spatially resolved spectra are shown for 
the galaxy center (blue), ±0.55 kpc (green) and ±1.1 kpc (or- 
ange), respectively. SDSS early- type galaxies described in Sec- 
tion 2.1 are plotted for comparison (grey dots). Dashed and 
solid curved-lines indicate empirical and theoretical demarcation 
lines between star-forming and AGN-host galaxies, respectively 
j Kauffmann et al.ll2003l : [Kewlev et al]|2001l '). Horizontal and verti- 
cal solid lines are demarcation lines for Seyfert and LINER galax- 
ies, indicati ng [Oiii] A5007/H/9 = 3 .0 and [Nil] A6584/Ho = 0.6, 
respectivelv l|Kauffmann et al.|[2003l ) . 

there seems no direct physical link between the emission- 
line gas and old stars. 

We note that the spatial profiles of the emission-line 
flux are asymmetric, while that of stellar continuum flux 
is symmetric. In Figure 6, for example, the line fluxes of 
[Oiii] A5007 and H/3 at —1.1 kpc aperture are a factor 
of two larger than those of the same emission lines at 
+ 1.1 kpc aperture (see also Figure 4). Such asymmetric 
spatial profiles of the emission lines further suggest that 
ionizing source may not be not physically linked to old 
stellar component. 

We investigate whether the distribution of young stel- 
lar population is asymmetric by measuring the 4000A 
break index D„(4000) as defined bv lBalogh et al.l (|1999f ) 
since £'„(4000) is an indicator of the presence of young 
stars. The D„(4000) indices measured from -1.1 kpc to 
1.1 kpc apertures are, respectively, [1.81, 1.86, 1.95, 1.87, 
1.80] with measurement error less than 0.01, showing 
that old stellar population is dominant for all aperture 
but with a frost of young stellar population. The sym- 
metric trend of the £'„(4000) indicates that the distribu- 
tion of young stellar population is also symmetric, further 
implying that ionizing source is not primarily linked to 
both old and young stellar populations. 

Figure 7 presents the AGN diagnostic diagram with 
flux ratios of [N ii] A6584 to Ha and [O ill] A5007 to H/3 
measured at each aperture. We flnd that both line ratios 
tend to decrease from the central to outer (±1.1 kpc) 
apertures, implying that ionizing power varies as a func- 
tion of the radial distance from the center. However, we 
note that the decrease of the [Nil] /Ha flux ratio can be 
interpreted as the decrease of gas density and/or metal- 




-1.1 -0.55 0.55 1.1 

Distance from the center (kpc) 

Fig. 8. — Kinematic properties of the nebular emission and stellar 
continuum along the slit. The upper and bottom panels show the 
spatial profiles of the line-of-sight velocity and velocity dispersion, 
respectively. Adopted color codes are the same as in Figure 6. 
The uncertainties of the line velocities are estimated by adding 
systematic and measurement errors in quadrature (upper panel) 
while the error bars of line dispersion indicate the measurement 
errors (lower panel). 

licity. To compare with the flux ratio measured from 
fiber observations, we plot the fiber-integrated emission- 
line ratios based on the SDSS spectrum. The integrated 
emission-line ratios are ^ 0.2 dex lower than those mea- 
sured at the center, implying that the emission-line ratios 
based on the fiber-integrated spectra can be significantly 
underestimated compared to the central emission-line ra- 
tios. 

4.2. Kinematics of Emission-line Gas 

Comparing the kinematic properties of emission-line 
gas and stellar component can provide an important clue 
on the nature of the ionizing sources. In Figure 8 we 
present the line-of-sight velocities and velocity disper- 
sions of each emission line and stellar component mea- 
sured at five different locations along the slit. The line- 
of-sight velocity measured from stellar lines at the galaxy 
center is 7218 ± 2 km s^^, which is consistent within the 
error with the SDSS redshift z = 0.0241 ± 0.0002, cor- 
responding to 7225 ± 60 km s~^. By adopting the line- 
of-sight velocity measured from stellar component at the 
galaxy center as a zero-point, we present the relative ve- 
locities of gas and stellar components. For emission lines. 
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we adopt the air wavelength provided by iMortod (|1991[ ) 
as a rest-frame wavelength to calculate the radial velocity 
of each emission line. 

We find several characteristics of the emission-line gas 
and stellar component. First, there is a rotating stel- 
lar component with a line-of-sight velocity of ^ ±90 
km s^^ at ^ ±1 kpc. Since the old stars are dominant 
in the luminosity-weighted galaxy spectra, the rotation 
is mainly represented by old stellar population. It is 
not unusual to find a rotating stellar component among 
early-type galaxies (e.g. EmscUcm ct al. 2007, 2011). For 
example, based on the integral field spectroscopic ob- 
servations of a sample of 48 nearby early- type galaxies, 
lEmsellem et al.l ()2007[ ) reported that all early-type galax- 
ies in their sample were either slow or fast rotator. 

Second, emission-line gas is counter-rotating with re- 
spect to stellar component, indicating that gas and stellar 
component are kinematically decoupled. The dynamical 
detachment implies that emission-line region may not be 
directly related to stellar component although we cannot 
rule out the stellar component as ionizing sources solely 
based on the dynamical decoupling. 

Third, [O ill] A5007 emission line is significantly 
blueshifted with respect to stellar component, by — 56 ± 
12 km at the center, while other emission lines do 
not show clear velocity offset. An artifact due to a 
misaligned slit position may cause the velocity offset 
([Comcrford ct al. 2009), however, the detected offset is 
not an artifact since the radial velocities of all other emis- 
sion lines except [O ill] A5007 are consistent with the ref- 
erence velocity of stellar lines. 

While the velocity offset of emission-line gas at the 
center of LINER galaxies has not been reportec0, it 
is rather common in Seyfert galaxies. Investigating 
[O II I] A5007 line for 65 loc al Seyfert galaxies, for exam- 
ple, iCrenshaw et al.l (|2010( ) summarized that 6% of the 
Seyferts display redshifted velocity offset with Av > +50 
km s~^, and 35% display blueshifted offset with Av < 
—50 km s~^. They provided a physical interpretation 
that the offset is due to the combined effect of gas out- 
flow of ionized gas in the NLR an d dust extinction in th e 
inner stellar disk (see Figure 7 of ICrenshaw etlll[201flh . 
Although this scenario was developed for Seyfert galax- 
ies, it is a viable solution for the velocity offset detected 
in the LINER galaxy, assuming that as a minor version of 
Seyfert galaxies, LINER galaxies also have mass outflow 
at the center (e.g., Kcwlcy ct al, 2006i) . Compared to 
other scenarios, e.g., infa lling small galaxy to the center 
of an AGN-host ga laxv ([Heckman et al.|[2009[ ) or align- 
ment of two AGNs (jShields et al.ll2009[ ) . the outflow sce- 
nario provides a simpler explanation for the velocity off- 
set of emission-line gas. Following the outflow scenario, 
it is probable that we only detected the [O ill] velocity 
offset because [O iii] ionization region is closer to the ac- 
cretion disk, hence, more affected by an outflow (e.g., 
iKomossa et"alll2008f ). 

Counter-rotating gas disks are often found in nearby 
early-type galaxies b ased o n longslit spectrosco py, e.g., 
NGC 4621 mended 11991 . NGC 4450 ^Rubin et al.l 
I1992D . and early- type galaxies of recent IFU surveys (e.g. 

^ Although IBarth et al.l II2008I ) reported that a radial velocity 
offset in one of LINER nuclei in NGC 3341, it is due to remote 
LINER nucleus from the center. 



lEmsellem et al.l 120071 120111 ). The observed velocity off- 
set of [Oili] may be explained by that outflow gas is 
mixed with the counter-rotating disk, producing the ve- 
locity offset with respect to stellar component and the 
rotation at the same time. It seems that spatial distri- 
bution and dynamics of NLR gas in the LINER galaxy is 
rather co mplex as siinilarly found in other local LINER 
galaxies (jPogge et al.ll2000l ). 

We investigate whether the LINER has an appropri- 
ate accretion rate to generate an outflow. Based on 
the empirical relationship between the mass of super- 
mass ive black hole (Mrh ) and stellar velocity disper- 
sion (iGiiltekin et al.l [2001 . we estimate black hole mass 
as « 1.1 X 10^ Mq from the stellar velocity dispersion 
178 ± 10 km s measured at the galaxy center. To 
measure bolometric luminosity (Lboi) we use the narrow 
Ha luminosity. Although it is more difficult to measure 
the Lboi of Type 2 AGNs than that of Type 1 AGNs 
due to the lack of broad emission-lines, the narrow Ha 
line can be used as a proxy for the broad emission-lines 
and the bolometric luminosity (Ho 2009). B y adoptin g 
the bolometric correction factor Cua 100 (iHol 120091) , 
we obtain L^oi ~ 6.0 x 10""^ erg s~^ from the narrow 
Ha line luminosity at the galaxy center. Combining the 
A/bh and L^oi estimates, we derive the Eddington ratio 
as ^ 4.2 X 10~^, which is simi lar to the Eddington ra- 
tios of nearby LINER galaxies ()Hcill2009() . In such a low 
Eddington ratio regime, AGN outflow is unlikely to be 
present due to the low radi ative pressure compar ed to 
the gravitational force fe.g.. iRaimundo et al.ll20ldl ). 

In contrast to the theoretical expectations, however, 
several observational studies reported LINER galaxies 
with a distinguishab le outflow (e.g., Walsh ct al. 200^; 
iMasegosa et al.ll2011[ ). For example, Walsh ct al. (200^ 
presented a clear evidence for nuclear outflow of several 
nearby LINER galaxies. In addition, NL R morpholo- 
gies o f nearby LINER galaxies are disturbed (jPogge et al.l 
l2000r ). implying that the disturbed morp hology is closely 
relate d to the outflow from the center (jMasegosa et al.l 
1201 Ih . Although it is still puzzling that how LINERs 
with low Eddington ratios can exhibit an outflow, the ob- 
served kinematics of the emission-line gas in the LINER 
galaxy seems to be connected with an outflow from the 
central AGN. 

5. DISCUSSIONS AND SUMMARY 

As a pilot study for investigating the nature of ioniz- 
ing sources in LINERs, we observed a LINER early- type 
galaxy, SDSS J091628.05-K420818.7 with Keck/LRIS to 
obtain spatially resolved spectra. By comparing with the 
PSF of the observations, we flnd that the emission-line 
region is slightly extend, but comparable to the seeing 
size, indicating that emission-line gas is centrally con- 
centrated and triggered by a central source. We find 
that emission-line gas is counter-rotating with respect to 
stellar component and a significant radial velocity offset 
of [O III] A5007 line. 

We suggest a scenario that gas outfiow in the NLR is 
ionized by a central source, presumably an active black 
hole, and mixed with gas in the counter-rotating disk. 
Though it is still difficult to disentangle the gas origin 
in this work, the gas is unlikely to be ionized and orig- 
inated from the evolved stars from the old stellar pop- 
ulation, since the gas is kinematically decoupled with 
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the old stellar pop ulation. The count er-rotating gas may 
have external (e.g. iRubin et al.lll992 ), or in ternal origins 
such as evolved stars fe.g. iHo et aLll2009D . and the gas 
may be causally linked with the existence of young stellar 
population. We cannot exclude the possibility that the 
young stars could be an ionizing source, but the contri- 
bution would be much weaker compared to that of AGN 
as the spatially concentrated profiles of emission- line flux 
indicate. 

Although some LINER galaxies may show an extended 
emission-line region, it is still possible that an active 
black hole is hidden in the integrated spectra. Thus, 
we cannot simply exclude the LINER galaxies identified 
in SDSS from the AGN population in the local universe. 
Detailed studies based on spatially-resolved spectra are 
necessary to reveal the true nature of the ionizing source 
of LINERs. 
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